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Abstract. Evidence of secular dynamical evolution for detached active bi- 
nary orbits are presented. First order decreasing rates of orbital angular mo- 
mentum (OAM), systemic mass (M = Mi + M 2 ) and orbital period of de- 
tached active binaries have been determined as J/J = —3.48 x 10 _10 yr -1 , 
M/M = -1.30 x KT^yr -1 and P/P = -3.96 x lCT^yr" 1 from the kinematical 
ages of 62 field detached systems. The ratio of d log J/d log M — 2.68 implies 
that either there are mechanisms which amplify AM loss <5 = 2.68 times with 
respect to isotropic AM loss of hypothetical isotropic winds or there exist exter- 
nal causes contributing AM loss in order to produce this mean rate of decrease 
for orbital periods. Various decreasing rates of OAM (d\ogJ/dt) and systemic 
mass (dlogM/dt) determine various speeds of dynamical evolutions towards a 
contact configuration. According to average dynamical evolution with S = 2.68, 
the fraction of 11, 23 and 39 per cent of current detached sample is expected to 
be contact system within 2, 4 and 6 Gyr respectively. 



Introduction 



Observational data and th eory of con tact binaries revised ex tensively bylMohnachil 
1 1981ft . IVilhnl (|!98lT l and iRncinskil (|l982T l. According to iBucinskil (|!98fiT l the 



most promising mechanism of formation for contact binaries involves the orbital 
angular momentum loss (AML) and the resulting orbita l decay of detached but 
close synchronized binaries. AML by magnetic brea king (ISchatzmaniri 959. Kraft 
119671 lMestel ll968^ became popular especially after ISkumanichl ((1972) who pre- 
sented observational evidence of decaying rotation rates for single stars. Mag- 
netic breaking and tidal locking have been considered as a main route to form 
W UMa-type contact binaries from the systems initially detached (cf. iHuand 
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Orbital period evolutions and time scale for forming contact systems from 
detached progenitors are predicted differently among the authors above. Since 
tidal locking is more effective at short periods, secular orbital period decreases 
were estimated slow at the beginning. Only the binaries of few days orbital 
periods were predicted to become contact systems within the order of ~ 10 9 
years. Shrinking orbits and related orbital period decreases, unfortunately, are 
not detectable on commonly used O-C diagrams formed by eclipse times due 
to: (1) short time-span covered by existing O-C data (at most 100 years), (2) 
large scattering of unevenly distributed O-C data, (3) existence of complicated 
larger amplitude shorter time-scale fluctuations by many different e ffects such 
as mass transfer, third companion an d/or magnetic cycles etc. (cf. iDemircanl 
I2000ll2002l : iKreiner. Kim &: Nhall200lh . Observed period decreases on O-C dia- 
grams, thus, cannot be counted as observational evidence for the secular orbital 
period decreases. 

There exists an opposing theory which predicts differen t scenario of con- 
tact binary formation by a fission process ([Roxburghl Il966h enabling contact 
binaries to form at the end of the pre-main sequence contraction. Consider- 
ing the brief lifetime of con tact stages (0.1 < t contact < 1 Gyr) estimated by 
iGuinan Bradstreetl (|l988h . this theory would fail to produce older population 
of W UMa binaries unless AM is conserved. However, conservation of AM would 
refute formation of contact systems from the detached progenitors. 

Debate on the formation mechanism c ontinues. Low space d ensity of contact 
binaries (0.2% in the solar neighborhood, iRucins kl 12001 12006T 1 implies limited 
lifetimes for contact stages and favors the formation mechanism from the de- 
tached progenitors. However, if tidal locking is efficient only at orbital periods 
comparable (~1 days) to contact systems, then the 0.2% space density become s 
too much that "contact binaries appear out of nowhere" (|Paczvnski et al\\2Q0(t ). 
Other mechanisms of OAM loss besides tidal locking must occur to account this 
observed space density in the solar neighborhood. 

This presentation aims to summarize the evidences of secular dynamical 
evolution for detached binary orbits rec ently obtained by u s from the kinemat- 
ics of chromospherically act ive binaries fearatas et al.ll2004T ) and from the kine- 
matics of W UMa systems ( Bilir et al.ll2005l ). Afte ra brief summary of or bital 



dynamics, mean dynamical evolution according to iDemircan et al. (2006) will 
be described on log J — log P, log M — log P and log J — log M diagrams. 



2. Evidence of secular decrease of orbital periods and masses 

First. IGuinan fe Bradstreetl (jl988f ) estimated a kinematical age of 8-10 Gyr for 
W UMa systems from their galactic space velocity dispersions. Theory of forma- 
tion from detached progenitors appears consiste nt even if t his age is compared 
to an earlier estimate of ~ 5 Gyr kinematical age (|Ekerll 19921 ) of possible progen- 
itors; the chromospherically active binaries (CAB). Being comparable to nuclear 
time-scale, the 8-10 Gyr age of W UMa systems implies that contac t syste ms 
either have been retaining their original AM according to iRoxburgbl (|l966l ) or 
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have been formed from detached binaries if orbits are unstable agai nst AM loss. 
Howev er, according to the space velocities and the dispersions by I Asian et al.l 
(1999), W UMa binaries are not older than RS CVn systems. 

Increased size of the samples (CAB &; W UMa) together with greatly im- 
prove d astrometric data (parallaxes k, proper motions) by Hipparcos ( Perrvman et al.l 
1997) motivated us to reanalyze the problem once more. First, identifying pos- 
sible members of young moving groups (MG) from 237 CAB binaries, we had 
spitted them into two sub samples of kinematically young (< 0.95 Gyr, N=95) 
and older field (3.86 Gyr, N=142) systems. 

MGs are kin ematically coherent groups of stars that share a common ori- 
gin. Egged (119941 ) defined a super cluster of stars gravitationally unbound in the 
solar neighborhood, but having a same kinematics while occupying the extended 
regions in the Galaxy. Therefore a MG, unlike well known open clusters covering 
only a limited s ky, can be obse rv ed at all dir ections. Kinematical criteria origi- 
nally defined by lEggenl (|l958al lbl Il989l . Il995h for determining possible m embers 
of the best-documented MGs are summarized bv iMontes et al.l (l2001al lbh. The 
basic idea is that a test stars space velocity vector must be equal and parallel, 
or at least with deviations smaller than the pre-determined limits, to the space 
velocity vector of a MG. The ages of MGs are known as open cluster ages from 
the main-sequence turn-off point. Ages of MG groups considered in this study 
are given in Table 1 together with kinematical parameters to identify them. 



Table 1. Parameters of best documented moving groups. 



Name 


Age 


(U, V, W) 


V T 


CP. 




(Myr) 


(km/s) 


(km/s) 


( Q \ n 


Local Association 


20 - 150 


(-11.6,-21.0,-11.4) 


26.5 


(5.98,-35.15) 


(Pleiades, a Per, M34, 






S Lyr, NGC 2516, IC2602) 










IC 2391 Superclustcr 


35 - 55 


(-20.6,-15.7,-9.1) 


27.4 


(5.82,-12.44) 


(IC 2391) 






Castor MG 


200 


(-10.7,-8.0,-9.7) 


16.5 


(4.75,-18.44) 


Ursa Major Group 


300 


(14.9,1.0,-10.7) 


18.4 


(20.55,-38.10) 


(Sirius Superclustcr) 










Hyadcs Superclustcr 


600 


(-39.7,-17.7,-2.4) 


43.5 


(6.40,6.50) 


(Hyadcs, Pracscpc) 











Dispersions of CAB and W UMa stars on the U — V diagrams are compared 
in Fig. 1, where possible MG members are also shown. U &V are space velocity 
components towards the Galactic center and in the direction of Galactic rotation. 
/.From the galactic space velocity dispersions, 3. 86 and 5.47 Gyr of kinematical 
ages were assigned to field CABs and W UMas ( Karatas et al.ll20"04l . iBilir et al.l 
2005). Although being smaller than earlier determinations, these ages too do 
not show a conflict with the theory of formation from detached progenitors. 
Moreover, the difference (1.61 G yr) between these ages was interpreted as a 
mean lifetime of contact stages bv lBilir et all (120051, . 

Formed by kinematical criteria, young groups allowed us to compare period 
and total mass histograms as well as other physical parameters between a young 
group and corresponding field systems in order to search further evidences of 
OAM loss. Since kinematical criteria determine only possible members for a 
MG, wrong identifications are always possible since the space velocity of a star 
coincidently may imply a membership. We believe number of wrong identifica- 
tions are small and thus unable to spoil statistics that evidences of dynamical 
evolution from detached to contact stages could still be found. 
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Figure 1. CAB (left), W UMa (right), Possible MG members (upper) are 
removed from the samples, then what left are called field systems (below). 
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Figure 2. Age dependent variations of OAM (J), period (P) and mass (M). 
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Being much younger than field systems, histograms of MG group CAB could 
be taken as initial distributions of field CAB. Co nsequently, it become s clear that 
smooth shaped initial distribution (see Fig. 6 of lKaratas et al.ll2004h is changed 
to a distribution with a peak at 10 days indicating the relative number of longer 
and shorter period systems were decreased. Secular OAM loss forcing orbital 
periods to decrease could explain the reduction of the relative number of systems 
with periods longer than 10 days. On the other hand, period decrease and radius 
increase as a result of dynamical and nuclear evolutions must have changed the 
missing short period systems into contact or semi-contact form thus they are no 
longer in the list of current CABs. We think the difference of period histograms 
between the MG group and the field CABs presents clear hints of OAM loss and 
secular decrease of orbital periods. 

It is not as clear, but similar kind of trends are also noticeable in the total 
mass (M = Mi + M 2 ) histograms (see Fig. 8 of lKaratas et al.ll2004^ . The tails 
rather than peaks support the prediction of mass decrease as a consequence of 
mass loss which carries away OAM of active binaries. The gradual decrease of the 
high mass tail of young group is changed to a sharper decrease in the older group. 
Similarly sharp decrease towards the less massive systems changed to rather a 
gradual decrease in the older population. Heterogeneity of samples containing 
giants (G), sub-giants (SG) and main-sequence (MS) systems and the evolu- 
tion into contact or semi-contact configurations complicat es the histograms and 
makes the interpretation of the peaks more difficult. See lKaratas et all 12004) 
for the details and further interpretations about sub-groups discriminating G, 
SG and MS systems. 

Increase of kinemat ical ages towards the short period systems (see Table 
5 of lKaratas et al.ll2004f ) could be taken as the clearest evidence of orbital pe- 
riod decrease as a consequence of OAM loss. Dynamical evolution has been 
further quantized with mean decreasing rates of OAM, total masses and peri- 
ods from the kinematical ages of carefully selec t ed de tached CABs with most 
reliable physical parameters bv iDemircan et all (|2006h . who preferred to form 
sub-samples according to OAM ranges. Mean OAM (J), period (P) and mass 
(M) of these sub groups versus the kinematical ages are displayed in Fig. 2. The 
decreases of J, P and M by age are obvious. As a first approximation linear 
lines were fitted by the least squares method and the inclinations were found to 
be -1.51 x 10~ 10 yr -1 , -1.72 x 10" 10 yr" 1 , and -5.65 x 10~ n yr" 1 . 

A linear change in the logarithmic scale implies a constant rate of change. 
Consequently, the mean relative decreasing rates J/J = —3.48 x 10 _10 yr _1 , 
P/P = -3.96 x lO^yr^and M/M = -1.30 x lO^yr" 1 were determined 
from kinematics to represent a mean dynamical evolution have been occurring 
among the detached CAB stars in the solar neighborhood. 

Similar histograms and tables were also produced for W UMa stars. How- 
ever, similar kind of interpretations are not pos sible for them because lifet ime 
in the cont act stage is too s hort, e.g. < 1 Gyr ( Guinan Sz Bradstreetlll988h or 
~ 1.6 Gyr (|Bilir et al.ll2005h . The MG group systems would disappear before 
reaching to an age of even the youngest group of field con tacts with a mean kine- 
matical age 3.21 Gyrs (See Table 6 of iBilir et al]l2005h . The four sub-groups 
formed from 97 field contact binaries according to orbital period ranges (ages 
3.21, 3.51, 7.14 and 8.89 Gyrs) dominate over 27 systems which are possible 
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MG members with ages less than 0.6 Gyrs. This fact implies that pre-contact 
dynamical evolution must exist dominantly among W UMa stars. Nevertheless, 
being out of the ordinary, di rect formation a t the end of pre-main sequence 
contraction by fission process (|Roxburerr] ll966) must also be occurring because 
ages of MG groups (< 0.6 Gyrs) give no space for pre-contact detached phases. 
Due to complexities of various lengths of pre-contact stages, it is not possible 
to evaluate the total mass and the period histograms for W UMa systems in a 
similar manner as CAB systems. 



3. Orbital dynamics 

The most basic definition of OAM (J) can be given as 



\M l + M 2 j \{l + q) 2 j y ' 

where / = (jf^j^) a 2 = ( ^ q y 2 )Ma 2 is moment of inertia and f2 = 2tt/P is 
angular speed for an orbital motion, thus J = 10. 

OAM (J) is needed dynamically to keep the orbital motion. Therefore, 
J and M are physical quantities which determine a unique period (P) and a 
unique size for the orbit as 



(l+ 9 ) 6 2vr J 3 (1+g) 4 J 2 m 

q3 G 2 M 5 ' a Gq 2 M3' { ) 

where the mass ratio (q = M2/M1 < 1) can be considered as an auxiliary 
parameter used in the definition of J. The size a = a± + a 2 represents the semi- 
major axis of a relative orbit of one star around the other. Stability of an orbit 
(dP = 0, da = 0) requires OAM and mass to be constant (dJ = 0, dM = 0) 
provided with no mass transfer (dq = 0). If there is no mass transfer, which 
must be true for detached binaries, it is obvious that OAM loss will cause an 
orbit to reduce its period and size. On the contrary, mass loss has an affect 
of increasing the period and size. Logarithmic derivatives of (2) give possible 
relative changes as 

dP _ 1-qdq dJ dM da _ 1 - q dq dJ dM 

~p " ~ 3 TT~ q J + 3 T " 5 m"' ~a ~ " 2 TT^7 + 2 T " 3 1T ' {3) 

Because M has higher power than J, the affect of mass loss would dominate. 
For example, in the case of same relative changes of OAM and mass (isotropic 
stellar winds, if dq = 0), dJ / J = dM/M according to (1), then 

dP _ dJ _ dM da _ dJ _ dM 

~P ~ ~ T ~~ ~W' ~a ~ ~T ~ ~~M' ^ ' 

which means mass loss and corresponding OAM loss will have a net effect on the 
orbit to increase both the period and the size. However, there could be additional 
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causes to increase relative OAM loss, e.g. 0AM loss of gravity waves, or stellar 
encounters in the galactic space, or a third body in an eccentric orbit around 
the binary system, or existence of an amplification mechanism as in some tidally 
locked binaries, in which tidal interactions transfer OAM to spinning components 
and AM is lost at the Alfven radius. After considering all possibilities, one has 
to compare the grant total relative OAM loss to the relative mass loss which 
could be expressed by a parameter S defined as 



which can be called dynamical parameter because dynamical respond of the 
orbit depends on the value of 5. Using this definition of 5, equation (3) becomes 



dP , 5 dJ dM da 3 dJ dM 

~P~ 5 ~ = ~ ~M~' T = (2 " 5 ] - = {26 ~ 3) 17' (6) 

if dq/q is neglected. Mass ratio change could be zero just because relative mass 
loss of components would be equal. Even if dq/q ^ 0, the term (1 — q)/(l + q) 
could be very small especially for high mass ratio (q ~ 1) systems. Thus, ignoring 
it in the first approximation is acceptable. 

For decreasing an orbital period, 5 > 5/3 is required. But, 5 > 3/2 is 
sufficient to shrink an orbit. If 3/2 < 5 < 5/3, orbital size decreases despite 
period is increasing. The size and the period of an orbit both increase if 5 < 3/2. 

Using the mean decreasing rates of OAM and mass from the kinematical 
ages of detached CAB stars, a mean value for the dynamical parameter (5) can 
be estimated for them as 

dJ dJ_ _o 4 o in-10 

A - -J- - Jdt - - 9 fiS (7\ 

— dM ~ dM — 1 on v in-10 ~ Z - UO y> 
-JT Mdt L - 6U X iU 

in the solar neighborhood. 



4. Mean dynamical evolution on diagrams 
4.1. The log J — logP diagram 

CAB and W UMa stars with available OAM are plotted on a log J — log P 
diagram (Fig. 3). The CAB stars containing giants, sub-giants, and main 
sequence and A & W type W UMa are indicated. W UMa stars are located 
at the lower left. Having larger masses and orbital periods, the CAB systems 
display a band like distribution on the right of W UMa stars as elongated from 
lower left to upper right. The constant total mass lines are computed with 
q = 0.88, which is the median value of mass ratios of CAB stars, using 

at which P is varied while chosen M values are fixed. 



8 



Eker, Demircan, Bilir, Karata§ 



V) 
CO 

o 



CO 

o 




**- * >. . ■» • 

4 * .. 



4M . 

»3m: 

2M ( 
1M, 



54 
53 
52 
51 
50 

-1 -0.5 0.5 1 1.5 2 2.5 3 3.5 
Log P (days) 

Figure 3. CAB stars containing giants (•), subgiants (+), main-sequence 
(*); W UMa stars A type (x), W type (o). Constant mass (M = Mi + M 2 ) 
lines ( ) are computed with q = 0.88 for CAB systems. 
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Figure 4. Upper boundary in (a) could be translated to (b) analytically. 
CAB (□), W UMa (o); MG (empty), field (filled). Young and old CABs are 
randomly mixed. Mean dynamical evolution (arrow), mean decreases (right 
sides of triangles) for 2 Gyrs. 
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Well defined smooth upper boundary of CAB systems appears as if tracing 
the path of dynamical evolution for them. OAM loss, mass loss and associated 
orbital period decrease would move a system from upper right to lower left as if 
parallel to the upper boundary. Finally, some systems would enter in the region 
of W UMa stars. However, this description is not quite correct according to 
the mean dynamical evolution with S = 2.68, which is marked by an arrow at 
the upper end of the upper boundary in Fig. 4b. The magnitude of the arrow 
symbolizes an interval of 2 Gyrs. The right sides of the triangle in the lower 
right show the amount of OAM loss and corresponding orbital period decrease. 

Fig. 4b plotted with symbols to indicate kinematically young (MG) and 
old (field) stars. The random distribution of the young and old CAB systems 
in implies that a detached system may start its dynamical evolution anywhere 
on the diagram in the region of CAB stars. So, it is not possible to distinguish 
younger and older systems according to their location on these diagrams unlike 
the nuclear evolution on H-R diagram which has a well defined starting point 
on the Zero Age Main Sequence (ZAMS). It is well known that the location of 
a star on H-R diagram indicates its stage of evolution as well as its age. 

4.2. The logM- logP diagram 

The log M — log P diagram resembles a similar type distribution as the log J — 
logP diagram. First, the upper boundary of CAB stars on the logM — logP 
diagram (Fig. 4a) was eye estimated and digitized by computer. Then using its 
numerical M, P points together with q = 0.88, the solid line on the log J— log P 
diagram was computed by (8). Since the computed line fits even better to the 
CAB upper boundary of log J— log P distribution, the upper boundaries of both 
diagrams in Fig. 4 are not independent and appears to be determined by the 
upper mass limits of orbital periods. Non-existence of CAB stars above those 
upper boundaries could be related to binary formation mechanism since initial 
Roche lobes may put limits on the masses of forming binaries. The binaries 
above, if they exist, they are not chromospherically active because being brighter 
they would have been easier to be noticed as CAB stars, then upper boundaries 
on those diagrams mark only the upper mass limit for a binary of a given period 
to have choromosphric activity. 

Moreover, If period decrease occurs because of OAM loss but without mass 
loss and mass transfer, a dynamical evolution would follow a path parallel to 
the constant mass lines in Fig. 3, which would carry binaries into the empty 
region. This requires for a CAB system to seize chromospheric activity if it 
moves into the empty region by such a process. Therefore, it is most likely the 
CAB upper boundary on log J — log P diagram represent a dynamical evolution 
with a minimum mass loss. Nevertheless, mean decreases determined from the 
kinematical ages of CAB stars shown by the right sides of the triangles in Fig. 4 
indicates that the direction of mean dynamical evolution neither is parallel nor 
towards the empty region as it is shown by the arrows at the upper ends of the 
upper boundaries in Fig. 4. 

4.3. The log J — logM diagram 

Because J and M are basic quantities to determine P & a and because OAM loss 
&l mass loss are parameters controlling the magnitude and direction of dynamical 
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Figure 5. Well defined borderline sharply separates detached and contact 
systems, (a) Symbols are like Fig. 3, (b) Symbols are like Fig. 4. Young and 
old CABs are randomly mixed. 



evolution, the log J — log M diagram is a natural choice to study dynamical 
evolution of orbits. Once, Fig. 5 is produced, a sharp separation between the 
detached and contact systems stroke to our attention. Goodness of separation is 
out striking that despite crowding along the border, there are only two systems 
(OO Aql, 5 Cap) on the wrong side, which could be due to a wrong identification 
of the state of being contact or just because of observational errors. Separation 
of such a degree does not occur on the diagrams discussed before. 

Marking several positions on the borderline between CABs and W UMas 
following quadratic equation was produced. 

log J = 0.522(log Mf + 1.664(log M) + 51.315, (9) 

where M is in solar units and J is in cgs. Physical significance of this line is 
that it marks the maximum OAM for a contact system to survive. It is like 
in single stars, spin AM has to be less than a certain value otherwise gravity 
cannot hold stellar mass together. It is same for contact binaries, if OAM is 
more than the value computed above, the contact configuration brakes. CAB 
systems (RT L ac, AR Mon, e UMi, R V Lib, BH CVn), which were eliminated 
from the list of lDemircan et ail (|2006T ) when the mean decreasing rates of J, M 
and P were determined (since they are filling or about to fill one of the Roche 
lobes that mass transfer possibly occurring in them), are marked on Fig. 6. 
Since those systems are not close to the contact border and scattered randomly 
all over in the detached region, the state of being detached and semi-detached 
must be same physically. 

A detached system must lose OAM to go into the region of contacts in 
order to be a contact binary. This possibility, however, very much depends on 
the position on the diagram together with the direction and the speed of the 
dynamical evolution, which could be very different from one system to another. 
The log J — log M diagram too gives no clue on individual dynamical evolutions 
since there is no information about initial positions. Random mix of young 



Dynamical Evolution of Detached Active Binary Orbits 



11 
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Figure 6. Loci of equal times to reach at the contact border. Mean dynam- 
ical evolution (hypothenuse) and mean decreases (right sides) for 2 Gyr. 



(MG) and old (field) systems (Fig. 5b) indicates that a same position could be 
belonging to both a young and an old system. 

Nevertheless, the contact border (eq. 9) could be helpful to determine con- 
tact binary candidates. The amount of mean losses (A J and AM) corresponding 
to 2, 4 and 6 Gyr are subtracted from the J and M values of the border. That 
is, the contact border is shifted accordingly. Dotted lines in Fig. 6 represent 
shifted borders that the systems between a dotted line and the contact border 
are the ones which are predicted to be contact systems within the time inter- 
vals indicated if their nuclear evolutions permit them to live untill the predicted 
times. After counting, it becomes clear that 11, 23 and 39 per cent of the current 
sample of CABs could pass over the contact border within the next 2, 4 and 6 
Gyr according to the mean dynamical evolution with 5 = 2.68. 

4.4. Period and size evolution of orbits 

It is possible to draw constant period lines using (8). On the other hand, 



J = . q , VGM 3 a (10) 
[1+qY 

can be used to compute constant orbital size lines similarly. From the statistics 
of present samples, median values q = 0.88 and q = 0.39 are found to represent 
CAB and W UMa stars. The constant orbital period and orbital size lines run 
almost parallel to the contact border and both PSza values decrease towards it 
(Fig. 7). Further decrease into the region of contacts is also clear. 

Constant period and constant size lines are sensitive to small q values. One 
can feel the sensitivity by comparing P = 1 day (or a = 3i? Q ) lines of q = 0.88 
and q = 0.39, also P = 0.2 days (or a = 1.2R & ) computed with q = 0.39 and 
q = 0.1. Decreasing periods (or sizes) towards the lower right is deceptive as 
if evolution to contact stage is occurring from upper left to lower right which 
is impossible since such an evolution requires a mass gain. Mass loss, however, 
changes the direction from vertically down (OAM loss only) towards to the lower 
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Figure 7. Period and size evolutions of orbits. Mean (hypotenuse) dynami- 
cal evolution and corresponding decreases (right sides of triangles). Evolution 
with 5 = 5/3 keeps P constant but 5 = 3/2 keeps a constant. 



left (if OAM and mass both are lost). Because the affect of mass loss dominates 
over OAM loss (eq. 3), there are lower limits; one for orbital periods and one for 
the orbital sizes. Both limits are indicated by the dotted lines in the triangles and 
the corresponding numerical values (d log J / d log M) in Fig. 7. Any dynamical 
evolution with a 8 smaller than those limits indicates an increase rather than a 
decrease on both P&a. 



5. Conclusions 

Since 5 = 2.68 for detached CAB in the solar neighbourhood is greater than both 
limits, (5/3 for P, 3/2 for a) orbits are shrinking and periods are decreasing. 
A well defined borderline sharply separating detached and contact systems are 
discovered empirically on log J — log M diagram. It is possible for a detached 
system to pass over the contact border and become a contact system by OAM 
loss. Not only P ~ 1 day period detached systems, but also some P ~ 10 days 
period systems can pass over this border according to the constant period lines in 
Fig. 7 and loci of equal times to reach at the contact border in Fig. 6. Significant 
number of current CABs (39 per cent) were found as canditates to pass over 
the contact border within the next 6 Gyrs according to their positions on the 
log J — log M diagram. Nevertheless, contact binary formation from detached 
progenitors is not the only mechanism to form contact systems. Although, it is 
rare, direct formation at the end of the pre-main sequence evolution must also 
be occurring. Statistical studies are encouraged to include these new findings if 
one can still say "contact b inaries appear out of nowhere" while they are already 
known to be rare systems ( Rucinskill2002l . I2006T ) . 
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